Passive scalar intermittency in low temperature helium flows 
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We report new measurements of turbulent mixing of tem- 
perature fluctuations in a low temperature helium gas exper- 
iment, spanning a range of microscale Reynolds number, R\, 
from 100 to 650. The exponents ^„ of the temperature struc- 
ture functions {\9{x + r) — 0{x)\") ~ r^" are shown to satu- 
rate to ^cxD — 1.45 ± 0.1 for the highest orders, n ~ 10. This 
saturation is a signature of statistics dominated by front-like 
structures, the cliffs. Statistics of the cliff characteristics are 
performed, particularly their width are shown to scale as the 
Kolmogorov length scale. 



The strong intermittency of a passive scalar field ad- 
vected by a turbulent flow has recently received con- 
siderable attention Two facets of this intermit- 
tency are the persistence of small scale anisotropy [|[^] 
and the anomalous scaling of the structure functions 
{\9{x + r) - e{x)Y') |,§. It is well established that the 
persistence of scalar gradient skewness arises from the 
ramp-and-cliff structures i.e high scalar jumps sep- 
arated by well mixed regions. The observed preferential 
alignment of cliffs with the large scale gradient ap- 
parently prevents a universal description of the odd order 
statistics, reflecting the asymetry of scalar fluctuations. 
However, the genericity of cliffs in "scalar turbulence" 
raises the issue of their influence on high-order scaling of 
even moments, and of the possible universality of the sus- 
pected saturation of high order exponents |9-11|. Labo- 
ratory experiments able to cover a wide range of Reynolds 
numbers in well controlled conditions appear to be crucial 
to address this point. We report in this Letter new mea- 
surements of turbulent mixing of a scalar fleld, namely 
temperature fluctuations, performed in a low tempera- 
ture helium gas experiment. Such measurements of tem- 
perature fluctuations, as a passively advected scalar field, 
are performed for the first time in low temperature he- 
lium gas , opening new and encouraging perspectives 
in investigation of turbulent mixing in high Reynolds 
number flows. 

The set-up we use is the same as the one described in 
Ref. p3| , with an additional apparatus to induce temper- 
ature fluctuations. The flow takes place in a cylindrical 
vessel and is driven by two rotating disks, as sketched in 
Fig. |l|. The disks, 20 cm in diameter and spaced 13.1 cm 
apart, are mounted with 6 radial blades. The cell is filled 
with helium gas, held at a controlled pressure. Tem- 
perature measurements are performed with a specially 
designed cold-wire thermometer in the thermal wake of 



a heated grid. The grid is made of a double nichrome 
wire net, stretched across a 4.0 x 4.5 cm frame. The wire 
diameter is 250 /Ltm, and the mesh size, M , is 2.0 mm. 
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FIG. 1. (a) Sketch of the experiment, (b) Upper view, in 
corotating mode. (1) Disks, (2) thermometer, (3) heated grid. 

The thermometer is located 20 mesh sizes downstream 
from the grid, at a distance of 2.2 cm from the wall. 
It has been designed as the velocity probes described in 
Ref. Ip^ . It consists in a thin carbon fiber, 7 /im in di- 
ameter, stretched across a 2 mm rigid frame. A metallic 
deposition of gold and silver of about 2000 A thick cov- 
ers the whole fiber, except for a 7 /im long central area. 
Temperature fluctuations are measured from resistance 
fluctuations, by the way of a constant current Wheat- 
stone bridge. The tuning of the current has received 
considerable attention, in order to ensure an high enough 
signal, without suffering from velocity contamination. In 
practice, the optimal value of the current is chosen in 
order to maximize the flatness factor of the temperature 
derivative. Temperature resolution is estimated, from 
the high-frequency white noise level, to around 100 fiK. 
The spatial resolution is limited by the probe size, which 
is at least twice smaller than the Kolmogorov scale in 
the present experiments. The frequency response is lim- 
ited by the constant current amplifler to 10 kHz, a value 
comfortably high to resolve the highest frequency of the 
temperature fluctuations. 

Since the same probe, with different electronic devices, 
is used as both thermometer and anemometer, we are 
able to perform temperature and velocity measurements 
at the same point in the same flow conditions. We have 
measured the various turbulence characteristics, such as 
the microscale Reynolds number Rx and the Kolmogorov 
scale 77, for the same disks rotation rates and kinematic 
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viscosity v, allowing us to characterize each temperature 
data set. We define here R\ = u'X/v the Reynolds num- 
ber based on the Taylor scale A = v! (y j e)^!'^ , where v! is 
the rms of the velocity fluctuations, e the mean dissipa- 
tion rate and r\ — (jr' jtf-l^. 

The disks rotation can be set in two modes, co- and 
counter-rotating, hereafter denoted by COR and CTR. 
In the corotating mode, the two disks are rotating in the 
same direction at the same speed, so that the flow can 
be thought as a solid-like rotation. The fluctuation rate 
u' /{U) is found to lie around 10 % in this value 
close to the one obtained in a previous experiment in a 
similar set-up In the counter- rotating mode, the 

turbulence intensity is much higher, but the speed ratio 
is tuned so that a strong mean advection remains at the 
height of the grid and the thermometer. This ensures 
a low enough fluctuation rate, between 13 and 25 %, 
allowing use of the Taylor hypothesis to convert temporal 
fluctuations into spatial ones. We have checked that the 
fluctuation rate is the same with and without the grid, 
meaning that the grid contribution to the turbulence can 
be considered as negligible. The only influence of the 
grid is injection of heated sheets in an already turbulent 
flow. Thus thermal turbulence originates in the mixing of 
these hot layers behind each wire, over a typical distance 
of M{U)/u' —1-2 cm downstream. 

The grid temperature has to be high enough to induce 
thermal fluctuations giving rise to an acceptable signal 
to noise ratio of at least 50 dB. The noise level of 100 /iK 
requires a rms signal of typically 6' = (0^)^/^ ~ 40 mK at 
the probe location, which is around 4 % of the grid over- 
heat. Since we are working with a closed flow, the heat 
transfer from the cell to the cryostat has to be balanced 
by the power supplied to the grid, in order to ensure ther- 
mal stationarity. This constraint limits the maximum du- 
ration of an experiment (up to 5 hours of continuous run), 
and the range of velocity and fluid density; in practice, 
values of R\ in the range 100-300 in the corotating case, 
and 200-650 in the counter-rotating case, are obtained. 
Although much lower than the highest R\ achievable in 
our set-up (up to 5000 ||l^ ) , wc hope that future noise im- 
provement of the constant current amplifier will increase 
the upper bound well beyond 650. 
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TABLE I. Typical experimental parameters, k = v/Pr is 
the thermal diffusivity, i/ the kinematic viscosity and Pr ~ 0.8. 
The sample size N* is expressed as N{U) / (27r/sr;) , where TV is 
the number of data points of the sample and fs the sampling 
rate. 



The temperature signal is filtered, then sampled and 
recorded on a 16 bit acquisition board (ITC-18 from In- 
struTech), at a sampling rate between 1 and 10 kHz, a 
value at least twice the low-pass filter frequency. Special 
attention has been paid to the low noise level of the whole 
amplifying and sampling channel. In order to achieve 
correct convergence of higher-order statistics, the sam- 
ple sizes are of the order of 10^ — 10^ Kolmogorov time 
scales long. Table | summarizes the characteristics of the 
typical data sets used here. 

In order to characterize the large scale of the tem- 
perature and velocity fluctuations, quantities of inter- 
est are the integral scales Ag and A^, deflned from the 
autocorrelation functions of temperature 9 and longitu- 
dinal velocity u. For our whole data set, we measure 
Ag = 7.0 ±1.5 mm and A„ = 9.6 ±1.0 mm, with 
no noticeable R\ dependence ||l^. This two values are 
rather close, indicating that kinetic energy and temper- 
ature variance are injected at the same large scale. The 
Prandtl number of the gas being close to unity, we are 
on the case where velocity and scalar fluctuations take 
place on the same range of scales. 
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FIG. 2. Probability density functions of the normalized 
temperature increments, for R\ — 650 (file #4 of Tab. ^. 
From the inner to the outer pdf, r/r] = 10*, 600, 30 and 3. 



The small scale intermittency of temperature fluctu- 
ations can be characterized through their distributions 
at different scales. Figure || shows the probability den- 
sity functions (pdf) of temperature increments A9{r) — 
6{x + r) — 9{x) for 4 different separations r, ranging from 
r/rj = 10*^ (large scale), 600, 30 (typical inertial scales) 
down to 3 (close to the dissipative scale). Each distri- 
bution has been normalized by its standard deviation 
a = (A6'(r)^)^/^. These pdf are strongly not self-similar, 
reflecting strong intermittency effects. The width of the 
tails is remarkable, showing scalar jumps of amplitude 
up to 40 times the standard deviation at the smallest 
scale. We can note the slight asymmetry of the distri- 
butions, linked to the well known property of small scale 
persistence of anisotropy . 
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FIG. 3. Temperature structure function exponents, 
normalized by ^2- O- = 280 in COR mode, 

A: Rx = 650 in CTR mode. The dashed line indicates the 
Corrsin-Obukhov scaling n/2. 
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FIG. 4. Ratio of two pdfs, at different inertial scales, 
ri — 45r; and r2 = 24677 {R\ = 280). The long dashed curve 
is a gaussian fit, and the horizontal dashed lines indicate the 
constant ratio of the pdf tails. Vertical arrows, at AO ~ 48' , 
indicate the dominant contribution to the 8th order moment. 

The evolution of these 

pdfs is characterized by the structure functions, defined 
as Sn{r) = {\d{x + r) — 9{x)\"), where angular brackets 
denote space (time) average. They are found to follow 
power law in terms of the scale r, namely Sn{r) ^ r^" . 
The measured structure function exponents divided 
by ^2, are shown in Fig. |, in the COR {Rx = 280) and 
CTR {Rx = 650) cases. The exponents are defined by 
plotting the compensated structure functions r~^"S'„(r) 
and tuning the value of the exponent to obtain a well 
defined plateau for inertial separations. This procedure 
allows to estimate the error bar for each order. Although 
the second order exponent presents some scatter (with 
a systematic increase from 0.45 to 0.65 with increasing 
R\), the normalization of the higher order exponents by 
^2 provides an excellent collapse for the different Rx ■ It is 
remarkable that, for comparable sample sizes, the highest 
available order is much lower at higher Reynolds num- 
ber: the fluctuations becomes much more intermittent 
and the convergence becomes poorer. For Rx — 650 we 
have to restrict to n < 8, whereas ?i = 12 can be achieved 



for Rx = 280. The two sets of exponents are consistent 
within error bars for n < 8, meaning that the large scale 
properties of the two flow configurations do not affect 
these inertial range statistics. 

The values of £,n/£,2 are found to strongly depart from 
the linear law n/2, and the gap increases with the order, 
a usual signature of inertial range intermittency of the 
passive scalar Furthermore, the exponents are found 
to increase extremely slowly with the order, strongly sug- 
gesting a saturation, for n around 10, at a value 



Coo = (2.3 ±0.1)^2 ^ 1.45 ±0.1 



(1) 



(where ^2 is considered equal to 2/3). This observation 
only relies on the COR data set, but we have noted that 
no deviation appears between the COR and CTR data 
for n < 8. A saturation of the structure function expo- 
nents is a signature of statistics dominated by shock-like 
structures, the thermal cliffs. A consequence of this sat- 
uration is a constant ratio of the far tails of the pdfs for 
inertial range separations. The ratio of two pdfs, for sep- 
arations ri and r2 well into the inertial range, is displayed 
in Fig. IJ. For temperature increments |A6'| > 46*', this 
ratio tends towards a constant. It must be noted that the 
jA^I ~ 49' parts of this distribution give the dominant 
contribution to the 8th order moment (see the vertical ar- 
rows), i.e they maximize the integrant |A0|®p(A0). This 
additional test confirms the observed saturation of high 
order exponents. 
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FIG. 5. Histograms of cliff width (a) and amplitude (b) for 
values of the threshold s = 3 (O), 7 (□) and 11 (o). The 
dashed curves on the figure (b) are log-normal fits. 

The observed trend towards a saturation of high or- 
der exponents motivates a detailed analysis of the cliffs. 
We have performed statistics of the cliff characteristics, 
defined as the strongest temperature gradients, singled 
out from time series of temperature fluctuations. We de- 
fine the cliffs from a simple threshold on the temperature 
derivative. 



\d0/dx\ > s {{de/dxY) 



2\l/2 



(2) 



where s is a non dimensional constant. Spatial deriva- 
tives are obtained from temporal ones using the Taylor 
hypothesis, and the gradients are estimated from finite 
difference over the smallest resolved separation. We de- 
fine the cliff amplitude 59 as the difference between the 
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two extrema surrounding the gradient satisfying (H) , and 
the chfF width A such that the temperature derivative 
takes values exceeding 90 % of its local maximum. His- 
tograms of amplitudes and widths are shown in Fig. 
for three values of the threshold s, at a Reynolds number 
Rx = 280. 

We first note that both amplitude and width his- 
tograms can be well fitted by log-normal distributions 
(given by a parabola in log- log coordinates). The mean 
amplitude of the cliffs is of order 6', and is found to 
increase with the threshold s. On the other hand, the 
cliff width remains constant for different thresholds. This 
means that the strength of the gradients depends mainly 
on the amplitude of the scalar jump, and not on its width. 
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of observed R\ (100-650), suggesting that the strongest 
cliffs remain concentrated over the smallest scale of the 
flow (see also Ref. ||l^). This observation may have im- 
portant consequences for processes such as reactive mix- 
ing or combustion, where the reaction rate is enhanced 
where concentration gradients are strong. Further in- 
sight into the cliffs contribution to the saturation of the 
high order structure function exponents needs a detailed 
study of their morphology and spatial distribution. Pre- 
liminary results Pq] , based on the waiting time between 
cliffs from the temperature time series, strongly suggests 
self-similar clustering for inertial separations. 
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and the European Commission's TMR programme. Con- 
tract no. ERBFMRXCT980175 "Intcrmittency" . 
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FIG. 6. Mean cliff width divided by the Kolmogorov length 
scale, as a function of R\. Q: COR, A: CTR. The dashed 
line is the average, (A)/r; ~ 13. 

A central issue concerning the cliffs is the Reynolds 
number dependence of their width. Figure |^ shows the 
mean cliff width (A), divided by the Kolmogorov scale 
r], for Rx ranging from 100 to 650. This plot extends a 
preliminary study performed only on the COR mode. 
We can see a well defined plateau, 

(A) = (13 ± 3)77. (3) 

The scatter is moderate, and probably reflects the diffi- 
culty to resolve properly the smallest scales. However, 
this plot confirms that the mean cliff width follows a 
R^ scaling. This law holds for values of the threshold 
s from 3 up to 40 in the case of the highest Rx (see Fig. |^) . 
We can note that the data from the two configurations 
overlap on the central region Rx = 200 — 300, suggesting 
that this small scale characteristics is not affected by the 
large scale properties of the flow. 

To summarize, we have performed for the first time 
measurements of turbulent mixing of temperature, con- 
sidered as a passive scalar field, in a low temperature he- 
lium experiment. The study of inertial range statistics of 
the temperature increments gives evidence of a saturation 
of the high order exponents, to a value ^00 — 1-45 ± 0.1. 
This observation reveals that inertial range statistics are 
dominated by the cliffs, concentrating large scalar jumps 
over small distance [p|-pT| . The cliff widths are shown to 
scale as the Kolmogorov length scale, in the whole range 
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